We are IntechOpen, 
the world's leading publisher of 


Open Access books 
Built by scientists, for scientists 


5.300 130,000 150M 


ailable International authors and editors Downloads 


Our author among the 


154 TOP 1% 12.2% 


Countries delivered to most cited s Contributors from top 500 universities 


Selection of our books indexed in the Book Citation Index 
in Web of Science™ Core Collection (BKCI) 


Interested in publishing with us? 
Contact book.department@intechopen.com 


Numbers displayed above are based on latest data collected. 
F information visit www.intechopen.com 


y 


18 


New Approaches for Theoretical Estimation of 
Mass Transfer Parameters in Both Gas-Liquid 
and Slurry Bubble Columns 


Stoyan NEDELTCHEV and Adrian SCHUMPE 
Institute of Technical Chemistry, TU Braunschweig 
Germany 


1. Introduction 


Bubble columns (with and without suspended solids) have been used widely as chemical 
reactors, bioreactors and equipment for waste water treatment. The key design parameters 
in bubble columns are: 

e gas holdup; 

e gas-liquid interfacial area; 

e volumetric liquid-phase mass transfer coefficient; 

e gas and liquid axial dispersion coefficients; 

Despite the large amount of studies devoted to hydrodynamics and mass transfer in bubble 
columns, these topics are still far from being exhausted. One of the essential reasons for 
hitherto unsuccessful modeling of hydrodynamics and mass transfer in bubble columns is 
the unfeasibility of a unified approach to different types of liquids. A diverse approach is 
thus advisable to different groups of gas-liquid systems according to the nature of liquid 
phase used (pure liquids, aqueous or non-aqueous solutions of organic or inorganic 
substances, non-Newtonian fluids and their solutions) and according to the extent of bubble 
coalescence in the respective classes of liquids. It is also necessary to distinguish consistently 
between the individual regimes of bubbling pertinent to a given gas-liquid system and to 
conditions of the reactor performance. 

The mechanism of mass transfer is quite complicated. Except for the standard air-water 
system, no hydrodynamic or mass transfer characteristics of bubble beds can be reliably 
predicted or correlated at the present time. Both the interfacial area a and the volumetric 
liquid-phase mass transfer coefficient kia are considered the most important design 
parameters and bubble columns exhibit improved values of these parameters (Wilkinson et 
al., 1992). For the design of a bubble column as a reactor, accurate data about bubble size 
distribution and hydrodynamics in bubble columns, mechanism of bubble coalescence and 
breakup as well as mass transfer from individual bubbles are necessary. Due to the complex 
nature of gas-liquid dispersion systems, the relations between the phenomena of bubble 
coalescence and breakup in bubble swarms and pertinent fundamental hydrodynamic 
parameters of bubble beds are still not thoroughly understood. 

The amount of gas transferred from bubbles into the liquid phase is determined by the 
magnitude of kra. This coefficient is an important parameter and its knowledge is essential 
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for the determination of the overall rate of chemical reaction in heterogeneous systems, i.e. 
for the evaluation of the effect of mass transport on the overall reaction rate. The rate of 
interfacial mass transfer depends primarily on the size of bubbles in the systems. The bubble 
size influences significantly the value of the mass transfer coefficient kı. It is worth noting 
that the effects of so-called tiny bubbles (d4«0.002 m) and large bubbles (d.20.002 m) are 
opposite. In the case of tiny bubbles, values of mass transfer coefficient increase rapidly as 
the bubble size increases. In the region of large bubbles, values of mass transfer coefficient 
decrease slightly as the bubble diameter increases. However, such conclusions have to be 
employed with caution. For the sake of correctness, it would therefore be necessary to 
distinguish strictly between categories of tiny and large bubbles with respect to the type of 
liquid phase used (e.g. pure liquids or solutions) and then to consider separately the values 
of liquid-phase mass transfer coefficient k, for tiny bubbles (with immobile interface), for 
large bubbles in pure liquids (mobile interface) and for large bubbles in solutions (limited 
interface mobility). 

The axial dispersion model has been extensively used for estimation of axial dispersion 
coefficients and for bubble column design. Some reliable correlations for the prediction of 
these parameters have been established in the case of pure liquids at atmospheric pressure. 
Yet, the estimations of the design parameters are rather difficult for bubble columns with 
liquid mixtures and aqueous solutions of surface active substances. 

Few sentences about the effect of high pressure should be mentioned. Hikita et al. (1980), 
Oztürk et al. (1987) and Idogawa et al. (1985a, b) in their gas holdup experiments at high 
pressure observed that gas holdup increases as the gas density increases. Wilkinson et al. 
(1994) have shown that gas holdup, kra and a increase with pressure. For design purposes, 
they have developed their own correlation which relates well kya and gas holdup. As the 
pressure increases, the gas holdup increases and the bubble size decreases which leads to 
higher interfacial area. Due to this reason, Wilkinson et al. (1992) argue that both a and kra 
will be underestimated by the published empirical equations. The authors suggest that the 
accurate estimation of both parameters requires experiments at high pressure. They 
proposed a procedure for estimation of these parameters on the basis of atmospheric results. 
It shows that the volumetric liquid-phase mass transfer coefficient increases with pressure 
regardless of the fact that a small decrease of the liquid-side mass transfer coefficient is 
expected. Calderbank and Moo-Young (1961) have shown that the liquid-side mass transfer 
coefficient decreases for smaller bubble size. The increase in interfacial area with increasing 
pressure depends partly on the relative extent to which the gas holdup increases with 
increasing pressure and partly on the decrease in bubble size with increasing pressure. 

The above-mentioned key parameters are affected pretty much by the bubble size 
distribution. In turn, it is controlled by both bubble coalescence and breakup which are 
affected by the physico-chemical properties of the solutions used. On the basis of dynamic 
gas disengagement experiments, Krishna et al. (1991) have confirmed that in the 
heterogeneous (churn-turbulent) flow regime a bimodal bubble size distribution exists: 
small bubbles of average size 5x10? m and fast rising large bubbles of size 5x10? m. 
Wilkinson et al. (1992) have proposed another set of correlations by using gas holdup data 
obtained at pressures between 0.1 and 2 MPa and extensive literature data. 

The flow patterns affect also the values of the above-mentioned parameters. Three different 
flow regimes are observed: 

e homogeneous (bubbly flow) regime; 

e transition regime; 
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e heterogeneous (churn-turbulent) regime. 

Under common working conditions of bubble bed reactors, bubbles pass through the bed in 
swarms. Kastanek et al. (1993) argue that the character of two-phase flow is strongly 
influenced by local values of the relative velocity between the dispersed and the continuous 
phase. On the basis of particle image velocimetry (PIV), Chen et al. (1994) observed three 
flow regimes: a dispersed bubble regime (homogeneous flow regime), vortical-spiral flow 
regime and turbulent (heterogeneous) flow regime. In the latter increased bubble-wake 
interactions are observed which cause increased bubble velocity. The vortical-spiral flow 
regime is observed at superficial gas velocity uc-0.021-0.049 m/s and is composed of four 
flow regions (the central plume region, the fast bubble flow region, the vortical-spiral flow 
region and the descending flow region) from the column axis to the column wall. 
According to Koide (1996) the vortical-spiral flow region might occur in the transition 
regime provided that the hole diameter of the gas distributor is small. Chen et al. (1994) 
have observed that in the fast bubble flow regime, clusters of bubbles or coalesced bubbles 
move upwards in a spiral manner with high velocity. The authors found that these bubble 
streams isolate the central plume region from direct mass exchange with the vortical-spiral 
flow region. In the heterogeneous flow regime, the liquid circulating flow is induced by 
uneven distribution of gas holdup. At low pressure in the churn-turbulent regime a much 
wider range of bubble sizes occurs as compared to high pressure. At low pressure there are 
large differences in rise velocity which lead to a large residence time distribution of these 
bubbles. In the churn-turbulent regime, frequent bubble collisions occur. 

Deckwer (1992) has proposed a graphical correlation of flow regimes with column diameter 
and uc. Another attempt has been made to determine the flow regime boundaries in bubble 
columns by using uc vs. gas holdup curve (Koide et al., 1984). The authors recommended 
that if the product of column diameter and hole size of the distributor is higher than 2x104 
m2, the flow regime is assumed to be a heterogeneous flow regime. In the bubble column 
with solid suspensions, solid particles tend to induce bubble coalescence, so the 
homogeneous regime is rarely observed. The transition regime or the heterogeneous regime 
is usually observed. 

In some works on mass transfer, the effects of turbulence induced by bubbles are 
considered. The flow patterns of liquid and bubbles are dynamic in nature. The time- 
averaged values of liquid velocity and gas holdup reveal that the liquid rises upwards and 
the gas holdup becomes larger in the center of the column. 

Wilkinson et al. (1992) concluded also that the flow regime transition is a function of gas 
density. The formation of large bubbles can be delayed to a higher value of superficial gas 
velocity (and gas holdup) when the coalescence rate is reduced by the addition of an 
electrolyte. Wilkinson and Van Dierendonck (1990) have demonstrated that a higher gas 
density increases the rate of bubble breakup especially for large bubbles. As a result, at high 
pressure mainly small bubbles occur in the homogeneous regime, until for very high gas 
holdup the transition to the churn-turbulent regime occurs because coalescence then 
becomes so important that larger bubbles are formed. The dependence of both gas holdup 
and the transition velocity in a bubble column on pressure can be attributed to the influence 
of gas density on bubble breakup. Wilkinson et al. (1992) argue that many (very) large 
bubbles occur especially in bubble columns with high-viscosity liquids. Due to the high rise 
velocity of the large bubbles, the gas holdup in viscous liquids is expected to be low, 
whereas the transition to the churn-turbulent regime (due to the formation of large bubbles) 
occurs at very low gas velocity. The value of surface tension also has a pronounced 


www.intechopen.com 


392 Mass Transfer in Multiphase Systems and its Applications 


influence on bubble breakup and thus gas holdup. When the surface tension is lower, fewer 
large bubbles occur because the surface tension forces oppose deformation and bubble 
breakup (Otake et al., 1977). Consequently, the occurrence of large bubbles is minimal due 
to bubble breakup especially in those liquids that are characterized with a low surface 
tension and a low liquid viscosity. As a result, relatively high gas holdup values are to be 
expected for such liquids, whereas the transition to the churn-turbulent regime due to the 
formation of large bubbles is delayed to relatively high gas holdup values. 


1.1 Estimation of bubble size 

The determination of the Sauter-mean bubble diameter ds is of primary importance as its 
value directly determines the magnitude of the specific interfacial area related to unit 
volume of the bed. All commonly recommended methods for bubble size measurement 
yield reliable results only in bubble beds with small porosity (gas holdupx0.06). The 
formation of small bubbles can be expected in units with porous plate or ejector type gas 
distributors. At these conditions, no bubble interference occurs. The distributions of bubble 
sizes yielded by different methods differ appreciably due to the different weight given to 
the occurrence of tiny bubbles. It is worth noting that the bubble formation at the orifice is 
governed only by the inertial forces. Under homogeneous bubbling conditions the bubble 
population in pure liquids is formed by isolated mutually non-interfering bubbles. 

The size of the bubbles leaving the gas distributor is not generally equal to the size of the 
bubbles in the bed. The difference depends on the extent of bubbles coalescence and break- 
up in the region above the gas distributor, on the distributor type and geometry, on the 
distance of the measuring point from the distributor and last but not least on the regime of 
bubbling. In coalescence promoting systems, the distribution of bubble sizes in the bed is 
influenced particularly by the large fraction of so-called equilibrium bubbles. The latter are 
formed in high porosity beds as a result of mutual interference of dynamic forces in the 
turbulent medium and surface tension forces, which can be characterized by the Weber 
number We. Above a certain critical value of We, the bubble becomes unstable and splits to 
bubbles of equilibrium size. On the other hand, if the primary bubbles formed by the 
distributor are smaller than the equilibrium size, they can reach in turbulent bubble beds the 
equilibrium size due to mutual collisions and subsequent coalescence. As a result, the mean 
diameter of bubbles in the bed again approaches the equilibrium value. In systems with 
suppressed coalescence, if the primary bubble has larger diameter than the equilibrium size, 
it can reach the equilibrium size due to the break-up process. If however the bubbles formed 
by the distributor are smaller than the equilibrium ones the average bubble size will remain 
smaller than the hypothetical equilibrium size as no coalescence occurs. Kastanek et al. 
(1993) argue that in the case of homogeneous regime the Sauter-mean bubble diameter d; 
increases with superficial gas velocity uc. 

The correct estimation of bubble size is a key step for predicting successfully the mass 
transfer coefficients. Bubble diameters have been measured by photographic method, 
electroresistivity method, optical-fiber method and the chemical-absorption method. 
Recently, Jiang et al. (1995) applied the PIV technique to obtain bubble properties such as 
size and shape in a bubble column operated at high pressures. 

In the homogeneous flow regime (where no bubble coalescence and breakup occur), bubble 
diameters can be estimated by the existing correlations for bubble diameters generated from 
perforated plates (Tadaki and Maeda, 1963; Koide et al., 1966; Miyahara and Hayashi, 1995) 
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or porous plates (Hayashi et al, 1975). Additional correlations for bubble size were 
developed by Hughmark (1967), Akita and Yoshida (1974) and Wilkinson et al. (1994). The 
latter developed their correlation based on data obtained by the photographic method in a 
bubble column operated between 0.1-1.5 MPa and with water and organic liquids. In 
electrolyte solutions, the bubble size is generally much smaller than in pure liquids 
(Wilkinson et al., 1992). 

In the transition regime and the heterogeneous flow regime (where bubble coalescence and 
breakup occur) the observed bubble diameters exhibit different values depending on the 
measuring methods. It is worth noting that the volume-surface mean diameter of bubbles 
measured near the column wall by the photographic method (Ueyama et al., 1980) agrees 
well with the predicted values from the correlation of Akita and Yoshida (1974). However, 
they are much smaller than those measured with the electroresistivity method and averaged 
over the cross-section by Ueyama et al. (1980). 

When a bubble column is operated at high pressures, the bubble breakup is accelerated due 
to increasing gas density (Wilkinson et al., 1990), and so bubble sizes decrease (Idogawa et 
al., 1985a, b; Wilkinson et al., 1994). Jiang et al. (1995) measured bubble sizes by the PIV 
technique in a bubble column operated at pressures up to 21 MPa and have shown that the 
bubble size decreases and the bubble size distribution narrows with increasing pressure. 
However, the pressure effect on the bubble size is not significant when the pressure is 
higher than 1.5 MPa. 

The addition of solid particles to liquid increases bubble coalescence and so bubble size. 
Fukuma et al. (1987a) measured bubble sizes and rising velocities using an electro-resistivity 
probe and showed that the mean bubble size becomes largest at a particle diameter of about 
0.2x103 m for an air-water system. The authors derived also a correlation. 

For pure, coalescence promoting liquids, Akita and Yoshida (1974) proposed an empirical 
relation for bubble size estimation based on experimental data from a bubble column 
equipped with perforated distributing plates. The authors argue that their equation is valid 
up to superficial gas velocities of 0.07 m/s. It is worth noting that Akita and Yoshida (1974) 
used a photographic method which is not very reliable at high gas velocities. The equation 
does not include the orifice diameter as an independent variable, albeit even in the 
homogeneous bubbling region this parameter cannot be neglected. 

For porous plates and coalescence suppressing media Koide and co-workers (1968) derived 
their own correlation. However, the application of this correlation requires exact knowledge 
of the distributor porosity. Such information can be obtained only for porous plates 
produced by special methods (e.g. electro-erosion), which are of little practical use, while 
they are not available for commonly used sintered-glass or metal plates. Kastanek et al. 
(1993) reported a significant effect of electrolyte addition on the decrease of bubble size. 
According to these authors, for the inviscid, coalescence-supporting liquids the ratio of 
Sauter-mean bubble diameter to the arithmetic mean bubble diameter is approximately 
constant (and equal to 1.07) within orifice Reynolds numbers in the range of 200-600. It is 
worth noting that above a certain viscosity value (higher than 2 mPa s) its further increase 
results in the simultaneous presence of both large and extremely small bubbles in the bed. 
Under such conditions the character of bubble bed corresponds to that observed for inviscid 
liquids under turbulent bubbling conditions. In such cases, only the Sauter-mean bubble 
diameter should be used for accurate bubble size characteristics. Kastanek et al. (1993) 
developed their own correlation (valid for orifice Reynolds numbers in between 200 and 
1000) for the prediction of Sauter-mean bubble diameter in coalescence-supporting systems. 
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According to it, the bubble size depends on the volumetric gas flow rate related to a single 
orifice, the surface tension and liquid viscosity. 

The addition of a surface active substance causes the decrease of Sauter-mean bubble diameter 
to a certain limiting value which then remains unchanged with further increase of the 
concentration of the surface active agent. It is frequently assumed that the addition of surface 
active agents causes damping of turbulence in the vicinity of the interface and suppression of 
the coalescence of mutually contacting bubbles. It is well-known fact that the Sauter-mean 
bubble diameters corresponding to individual coalescent systems differ only slightly under 
turbulent bubbling conditions and can be approximated by the interval 6-7x103 m. 


1.2 Estimation of gas holdup 

Gas holdup is usually expressed as a ratio of gas volume Vc to the overall volume (Vc* Vi). 
It is one of the most important parameters characterizing bubble bed hydrodynamics. The 
value of gas holdup determines the fraction of gas in the bubble bed and thus the residence 
time of phases in the bed. In combination with the bubble size distribution, the gas holdup 
values determine the extent of interfacial area and thus the rate of interfacial mass transfer. 
Under high gas flow rate, gas holdup is strongly inhomogeneous near the gas distributor 
(Kiambi et al., 2001). 

Gas holdup correlations in the homogeneous flow regime have been proposed by Marrucci 
(1965) and Koide et al. (1966). The latter is applicable to both homogeneous and transition 
regimes. It is worth noting that the predictions of both equations agree with each other very 
well. Correlations for gas holdup in the transition regime are proposed by Koide et al. (1984) 
and Tsuchiya and Nakanishi (1992). Hughmark (1967), Akita and Yoshida (1973) and Hikita 
et al. (1980) derived gas holdup correlations for the heterogeneous flow regime. The effects 
of alcohols on gas holdup were discussed and the correlations for gas holdups were 
obtained by Akita (1987a) and Salvacion et al. (1995). Koide et al. (1984) argues that the 
addition of inorganic electrolyte to water increases the gas holdup by 20-30 % in a bubble 
column with a perforated plate as a gas distributor. Akita (1987a) has reported that no 
increase in gas holdup is recognized when a perforated plate of similar performance to that 
of a single nozzle is used. Oztürk et al. (1987) measured gas holdups in various organic 
liquids in a bubble column, and have reported that gas holdup data except those for mixed 
liquids with frothing ability are described well by the correlations of Akita and Yoshida 
(1973) and Hikita et al. (1980). Schumpe and Deckwer (1987) proposed correlations for both 
heterogeneous flow regime and slug flow regime in viscous media including non- 
Newtonian liquids. Addition of a surface active substance (such as alcohol) to water inhibits 
bubble coalescence and results in an increase of gas holdup. Grund et al. (1992) applied the 
gas disengagement technique for measuring the gas holdup of both small and large bubble 
classes. Tap water and organic liquids were used. The authors have shown that the 
contribution of small class bubbles to kra is very large, e.g. about 68 % at uc=0.15 m/s in an 
air-water system. Grund et al. (1992) suggested that a rigorous reactor model should 
consider two bubble classes with different degrees of depletion of transport component in 
the gas phase. Muller and Davidson (1992) have shown that small-class bubbles contribute 
20-50 % of the gas-liquid mass transfer in a column with highly viscous liquid. Addition of 
solid particles to liquid in a bubble column reduces the gas holdup and correlations of gas 
holdup valid for transition and heterogeneous flow regimes were proposed by Koide et al. 
(1984), Sauer and Hempel (1987) and Salvacion et al. (1995). 
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Wilkinson et al. (1992) have summarized some of the most important gas holdup 
correlations and have discussed the role of gas density. The authors reported also that at 
high pressure gas holdup is higher (especially for liquids of low viscosity) while the average 
bubble size is smaller. Wilkinson et al. (1992) determined the influence of column 
dimensions on gas holdup. Kastanek et al. (1993) reported that at atmospheric pressure the 
gas holdup is virtually independent of the column diameter provided that its value is larger 
than 0.15 m. This information is critical to scale-up because it determines the minimum scale 
at which pilot-plant experiments can be implemented to estimate the gas holdup (and mass 
transfer) in a large industrial bubble column. Wilkinson et al. (1992) reached this conclusion 
for both low and high pressures and in different liquids. 

Wilkinson et al. (1992) argues that the gas holdup in a bubble column is usually not uniform. 
In general, three regions of different gas holdup are recognized. At the top of the column, 
there is often foam structure with a relatively high gas holdup, while the gas holdup near 
the sparger is sometimes measured to be higher (for porous plate spargers) and sometimes 
lower (for single-nozzle spargers) than in the main central part of the column. The authors 
argue that if the bubble column is very high, then the gas holdup near the sparger and in the 
foam region at the top of the column has little influence on the overall gas holdup, while the 
influence can be significant for low bubble columns. The column height can influence the 
value of the gas holdup due to the fact that liquid circulation patterns (that tend to decrease 
the gas holdup) are not fully developed in short bubble columns (bed aspect ratio<3). All 
mentioned factors tend to cause a decrease in gas holdup with increasing column height. 
Kastanek et al. (1993) argues that this influence is negligible for column heights greater than 
1-3 m and with height to diameter ratios above 5. 

Wilkinson et al. (1992) have shown that the influence of the sparger design on gas holdup is 
negligible (at various pressures) provided the sparger hole diameters are larger than 
approximately 1-2x103 m (and there is no maldistribution at the sparger). In high bubble 
columns, the influence of sparger usually diminishes due to the ongoing process of bubble 
coalescence. Wilkinson et al. (1992) argue that the relatively high gas holdup and mass 
transfer rate that can occur in small bubble columns as a result of the use of small sparger 
holes will not occur as noticeably in a high bubble column. In other words, a scale-up 
procedure, in which the gas holdup, the volumetric mass transfer coefficient and the 
interfacial area are estimated on the basis of experimental data obtained in a pilot-plant 
bubble column with small dimensions (bed aspect ratio<5, D.<0.15 m) or with porous plate 
spargers, will in general lead to a considerable overestimation of these parameters. Shah et 
al. (1982) reported many gas holdup correlations developed on the basis of atmospheric data 
and they do not incorporate any influence of gas density. 

In the case of liquid mixtures, Bach and Pilhofer (1978), Godbole et al. (1982) and Khare and 
Joshi (1990) determined that gas holdup does not decrease if the viscosity of water is 
increased by adding glycerol, carboxymethyl cellulose (CMC) or glucose but passes through 
a maximum. Wilkinson et al. (1992) assumes that this initial increase in gas holdup is due to 
the fact that the coalescence rate in mixtures is lower than in pure liquids. The addition of an 
electrolyte to water is known to hinder coalescence with the result that smaller bubbles 
occur and a higher gas holdup than pure water. 

The addition of solids to a bubble column will in general lead to a small decrease in gas 
holdup (Reilly et al., 1986) and the formation of larger bubbles. The significant increase in 
gas holdup that occurs in two-phase bubble columns (due to the higher gas density) will 
also occur in three-phase bubble columns. A temperature increase leads to a higher gas 
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holdup (Bach and Pilhofer, 1978). A change in temperature can have an influence on gas 
holdup for a number of reasons: due to the influence of temperature on the physical 
properties of the liquid, as well as the influence of temperature on the vapor pressure. 

Akita and Yoshida (1973) proposed their own correlation for gas holdup estimation. The 
correlation can be safely employed only within the set of systems used in the author's 
experiments, i.e. for systems air (O2, He, CO2)-water, air-methanol and air-aqueous solutions 
of glycerol. The experiments were carried out in a bubble column 0.6 m in diameter. The 
clear liquid height ranged between 1.26 and 3.5 m. It is worth noting that the effect of 
column diameter was not verified. Hikita and co-workers (1981) proposed another complex 
empirical relation for gas holdup estimation based on experimental data obtained in a small 
laboratory column (column diameter-0.1 m, clear liquid height=0.65 m). Large set of gas- 
liquid systems including air-(H»?, CO», CH, C3Hs, No)-water, as well as air-aqueous 
solutions of organic liquids and electrolytes were used. For systems containing pure organic 
liquids the empirical equation of Bach and Pilhofer (1978) is recommended. The authors 
performed measurements in the systems air-alcohols and air-halogenated hydrocarbons 
carried out in laboratory units 0.1-0.15 m in diameter, at clear liquid height > 12 m. 
Hammer and co-workers (1984) proposed an empirical correlation valid for pure organic 
liquids at low superficial gas velocities. The authors pointed out that there is no any relation 
in the literature that can express the dependence of gas holdup on the concentration in 
binary mixtures of organic liquids. The effect of the gas distributor on gas holdup can be 
important particularly in systems with suppressed bubble coalescence. The majority of 
relations can be employed only for perforated plate distributors, while considerable increase 
of gas holdup in coalescence suppressing systems is observed in units with porous 
distributors. Kastanek et al. (1993) argue that the distributor geometry can influence gas 
holdup in turbulent bubble beds even in coalescence promoting systems at low values of 
bed aspect ratio and plate holes diameter. 

The gas holdup increases with decreasing surface tension due to the lower rise velocity of 
bubbles. The effect of surface tension in systems containing pure liquids is however only 
slight. Gas holdup is strongly influenced by the liquid phase viscosity. However, the effect 
of this property is rather controversial. The effect of gas phase properties on gas holdup is 
generally of minor importance and only gas viscosity is usually considered as an important 
parameter. Large bubble formation leads to a decrease in the gas holdup. Kawase et al. 
(1987) developed a theoretical correlation for gas holdup estimation. Godbole et al. (1984) 
proposed a correlation for gas holdup prediction in CMC solutions. 

Most of the works in bubble columns dealing with gas holdup measurement and prediction 
are based on deep bubble beds (Hughmark, 1967; Akita and Yoshida, 1973; Kumar et al., 
1976; Hikita et al., 1980; Kelkar et al., 1983; Behkish et al., 2007). A unique work concerned 
with gas holdup £c under homogeneous bubbling conditions was published by Hammer et 
al. (1984). The authors presented an empirical relation valid for pure organic liquids at 
uc<0.02 m-s. Idogawa et al. (1987) proposed an empirical correlation for gas densities up to 
121 kg-m? and uc values up to 0.05 m.s. Kulkarni et al. (1987) derived a relation to 
compute zc in the homogeneous flow regime in the presence of surface-active agents. By 
using a large experimental data set, Syeda et al. (2002) have developed a semi-empirical 
correlation for gc prediction in both pure liquids and binary mixtures. PoSarac and Tekić 
(1987) proposed a reliable empirical correlation which enables the estimation of gas holdup 
in bubble columns operated with dilute alcohol solutions. A number of gas holdup 
correlations were summarized by Hikita et al. (1980). Recently, Gandhi et al. (2007) have 
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proposed a support vector regression-based correlation for prediction of overall gas holdup 
in bubble columns. As many as 1810 experimental gas holdups measured in various 
gas-liquid systems were satisfactorily predicted (average absolute relative error: 12.1%). The 
method is entirely empirical. 

In the empirical correlations, different dependencies on the physicochemical properties and 
operating conditions are implicit. This is primarily because of the limited number of liquids 
studied and different combinations of dimensionless groups used. For example, the gas 
holdup correlation proposed by Akita and Yoshida (1973) can be safely employed only 
within the set of systems used in the authors' experiments (water, methanol and glycerol 
solutions). The effect of column diameter D. was not verified and the presence of this 
parameter in the dimensionless groups is thus only formal. In general, empirical correlations 
can describe sg data only within limited ranges of system properties and working 
conditions. In this work a new semi-theoretical approach for gc prediction is suggested 
which is expected to be more generally valid. 


1.3 Estimation of volumetric liquid-phase mass transfer coefficient 

The volumetric liquid-phase mass transfer coefficient is dependent on a number of variables 
including the superficial gas velocity, the liquid phase properties and the bubble size 
distribution. The relation for estimation of kta proposed by Akita and Yoshida (1974) has 
been usually recommended for a conservative estimate of kra data in units with perforated- 
plate distributors. The equation of Hikita and co-workers (1981) can be alternatively 
employed for both electrolytes and non-electrolytes. However, the reactor diameter was not 
considered in their relation. Hikita et al. (1981), Hammer et al. (1984) and Merchuk and Ben- 
Zvi (1992) developed also a correlation for prediction of the volumetric liquid-phase mass 
transfer coefficient kra. 

Calderbank (1967) reported that values of kj, decrease with increasing apparent viscosity 
corresponding to the decrease in the bubble rise velocity which prolongs the exposure time 
of liquid elements at the bubble surface. The ki, value for the frontal area of the bubble is 
higher than the one predicted by the penetration theory and valid for rigid spherical bubbles 
in potential flow. The rate of mass transfer per unit area at the rear surface of spherical-cap 
bubbles in water is of the same order as over their frontal areas. For more viscous liquids, 
the equation from the penetration theory gives higher values of k, than the average values 
observed over the whole bubble surface which suggests that the transfer rate per unit area at 
the rear of the bubble is less than at its front. 

Calderbank (1967) reported that the increase of the pseudoplastic viscosity reduces the rate 
of mass transfer generally, this effect being most substantial for small bubbles and the rear 
surfaces of large bubbles. The shape of the rear surface of bubbles is also profoundly 
affected. Evidently these phenomena are associated with the structure of the bubble wake. 
In the case of coalescence promoting liquids, almost no differences have been reported 
between kra values determined in systems with large- or small-size bubble population. For 
coalescence suppressing systems, it is necessary to distinguish between aqueous solutions of 
inorganic salts and aqueous solutions of surface active substances in which substantial 
decrease of surface tension occurs. Values of kya reported in the literature for solutions of 
inorganic salts under conditions of suppressed bubble coalescence are in general several 
times higher than those for coalescent systems. On the other hand, kra values observed in 
the presence of surface active agents can be higher or lower than those corresponding to 
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pure water. No quantitative relations are at present available for prediction of kra in 
solutions containing small bubbles. The relation of Calderbank and Moo-Young (1961) is 
considered the best available for the prediction of kų values. It is valid for bubble sizes 
greater than 2.5x10? m and systems water-oxygen, water-CO2 and aqueous solutions of 
glycol or polyacrylamide-CO». For small bubbles of size less than 2.5x103 m in systems of 
aqueous solutions of glycol-CO», aqueous solutions of electrolytes-air, waxes-H» these 
authors proposed another correlation. An exhaustive survey of published correlations for 
kia and kr was presented by Shah and coworkers (1982). The authors stressed the important 
effect of both liquid viscosity and surface tension. Kawase and Moo-Young (1986) proposed 
also an empirical correlation for kta prediction. The correlation developed by Nakanoh and 
Yoshida (1980) is valid for shear-thinning fluids. 

In many cases of gas-liquid mass transfer in bubble columns, the liquid-phase resistance to 
the mass transfer is larger than the gas-phase one. Both the gas holdup and the volumetric 
liquid-phase mass transfer coefficient kia increase with gas velocity. The correlations of 
Hughmark (1967), Akita and Yoshida (1973) and Hikita et al. (1981) predict well kra values 
in bubble columns of diameter up to 5.5 m. Öztürk et al. (1987) also proposed correlation for 
kra prediction in various organic liquids. Suh et al. (1991) investigated the effects of liquid 
viscosity, pseudoplasticity and viscoelasticity on kia in a bubble column and they 
developed their own correlation. In highly viscous liquids, the rate of bubble coalescence is 
accelerated and so the values of kia decrease. Akita (1987a) measured the kra values in 
inorganic aqueous solutions and derived their own correlation. Addition of surface-active 
substances such as alcohols to water increases the gas holdup, however, values of ka in 
aqueous solutions of alcohols become larger or smaller than those in water according to the 
kind and concentration of the alcohol (Salvacion et al., 1995). Akita (1987b) and Salvacion et 
al. (1995) proposed correlations for kya prediction in alcohol solutions. 

The addition of solid particles (with particle size larger than 10 um) increases bubble 
coalescence and bubble size and hence decreases both gas holdup and kra. For these cases, 
Koide et al. (1984) and Yasunishi et al. (1986) proposed correlations for kta prediction. Sauer 
and Hempel (1987) proposed kia correlations for bubble columns with suspended particles. 
Sada et al. (1986) and Schumpe et al. (1987) proposed correlations for kia prediction in 
bubble columns with solid particles of diameter less than 10 um. Sun and Furusaki (1989) 
proposed a method to estimate kra when gel particles are used. Sun and Furusaki (1989) and 
Salvacion et al. (1995) showed that kra decreases with increasing solid concentration in gel- 
particle suspended bubble columns. Salvacion et al. (1995) showed that the addition of 
alcohol to water increases or decreases kra depending on the kind and concentration of the 
alcohol added to the water and proposed a correlation for kia including a parameter of 
retardation of surface flow on bubbles by the alcohol. 


1.4 Estimation of liquid-phase mass transfer coefficient 

The liquid-phase mass transfer coefficients kų are obtained either by measuring kia, gas 
holdup and bubble size or by measuring kra and a with the chemical absorption method. 
Due to the difficulty in measuring distribution and the averaged value of bubble diameters 
in a bubble column, predicted values of ky by existing correlations differ. Hughmark (1967), 
Akita and Yoshida (1974) and Fukuma et al. (1987b) developed correlations for k 
prediction. In the case of slurry bubble columns, Fukuma et al. (1987b) have shown that the 
degrees of dependence of kų on both bubble size and the liquid viscosity are larger than 
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those in a bubble column. Schumpe et al. (1987) have shown that low concentrations of high 
density solids of size less than 10 um increase kų by a hydrodynamic effect on the liquid film 
around the bubbles. 

For pure liquids and large bubbles (d.20.002 m), Higbie's (1935) relation based on the 
penetration theory of mass transfer can be used as the first approximation yielding 
qualitative information on the effect of fundamental physico-chemical parameters (viscosity, 
density, surface tension) on k values. All these parameters influence both the size of the 
bubbles (and consequently also their ascending velocity) and the hydrodynamic situation at 
the interface (represented by an appropriate value of liquid molecular diffusivity). Kastanek 
et al. (1993) proposed their own correlation for calculation of k;. 

Values of kų decrease with increasing apparent viscosity corresponding to the decrease in 
bubble rise velocity which prolongs the exposure time of liquid elements at the bubble 
surface. According to Calderbank (1967), k; for the frontal area is 1.13 times higher than the 
one predicted by the penetration theory and valid for spherical bubbles in potential flow. In 
the case of water, the rate of mass transfer per unit area at the rear surface of spherical-cap 
bubbles is of the same order as over their frontal areas. For more viscous liquids, the transfer 
rate per unit area at the rear of the bubble is less than at its front. 

Calderbank (1967) reported that in general the increase of pseudoplastic viscosity reduces 
the rate of mass transfer, this effect being most substantial for small bubbles and the rear 
surfaces of large bubbles. The shape of the rear surface of bubbles is also profoundly 
affected. According to Calderbank (1967), these phenomena are associated with the 
structure of the bubble wake. Calderbank and Patra (1966) have shown experimentally that 
the average ky obtained during the rapid formation of a bubble at a submerged orifice is less 
than the value observed during its subsequent ascent. According to the authors, this is a 
consequence of the fact that if the rising bubbles are not in contact with each other the mean 
exposure time of liquid elements moving round the surface of a rising bubble must be less 
than the corresponding exposure time during its formation. 

Large bubbles (d, >2.5x103 m) have greater mass transfer coefficients than small bubbles 
(d,«2.5x103 m). Small "rigid sphere" bubbles experience friction drag, causing hindered 
flow in the boundary layer sense. Under these circumstances the mass transfer coefficient is 
proportional to the two-thirds power of the diffusion coefficient (Calderbank, 1967). For 
large bubbles (22.5x103 m) form drag predominates and the conditions of unhindered flow 
envisaged by Higbie (1935) are realized. The author assumed unhindered flow of liquid 
round the bubble and destruction of concentration gradients in the wake of the bubble. 
Griffith (1960) suggested that the mass transfer coefficient for the region outside a bubble 
may be computed if one knows the average concentration of solute in the liquid outside the 
bubble, the solute concentration at the interface and the rate of solute transfer. Leonard and 
Houghton (1963) reported that the ky values for pure carbon dioxide bubbles dissolving in 
water is proportional to the square of the instantaneous bubble radius for diameters in the 
range 6-11x10? m where the rise velocity appeared to be independent of size. Leonard and 
Houghton (1961) found that for bubbles with diameters below 6x10? m mass transfer seems 
to have an appreciable effect upon the velocity of rise, indicating that surface effects 
predominate in this range of sizes. Hammerton and Garner (1954) argue that there is a 
simple hydrodynamic correspondence between bubble velocity and mass transfer rate. 
According to Leonard and Houghton (1963) kı is not only a function of bubble diameter but 
is also a function of the distance from the point of release. The variation of kų with distance 
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from the release point indicates that the rate is a function of time after release or some other 
related variable such as bubble size or hydrostatic pressure. Baird and Davidson (1962) 
observed a time dependence for carbon dioxide bubbles in water, but only for bubbles larger 
than 25x10- m in diameter, the explanation being that the time dependence was due to the 
unsteady state eddy diffusion into the turbulent wake at the rear of the bubble. Davies and 
Taylor (1950) developed a relation for k; prediction in potential flow around a spherical-cap 
bubble. The authors argue that the bubble shape becomes oblate spheroidal for bubble sizes 
below 15x102 m. 

Leonard and Houghton (1963) reported that the effect of inert gas is to reduce somewhat 
the mass transfer rate by about 20-40 % and to introduce more scatter in the calculated 
values of ki, presumably because of the smaller volume changes. Gas circulation is also 
involved. The effect of an inert gas is to reduce the specific absorption rate, presumably by 
providing a gas-film resistance that may be affected by internal circulation. 

Leonard and Houghton (1963) argue that there is a detectable decrease of kų with increasing 
distance from the release point during absorption, the reverse appearing to be true for 
desorption. The addition of surfactant can reduce mass transfer without affecting the rise 
velocity. Mass transfer from single rising bubbles is governed to a large extent by surface 
effects, particularly at the smaller sizes. 

The theory of isotropic turbulence can be used also for ky prediction (Deckwer, 1980). The 
condition of local isotropy is frequently encountered. The theory of local isotropy gives 
information on the turbulent intensity in the small volume around the bubble. Turbulent 
flow produces primary eddies which have a wavelength or scale of similar magnitude to the 
dimensions of the main flow stream. These large primary eddies are unstable and 
disintegrate into smaller bubbles until all their energy is dissipated by viscous flow. When 
the Reynolds number of the main flow is high most of the kinetic energy is contained in the 
large eddies but nearly all of the dissipation occurs in the smallest eddies. If the scale of the 
main flow is large in comparison with that of the energy-dissipating eddies a wide spectrum 
of intermediate eddies exist which contain and dissipate little of the total energy. The large 
eddies transfer energy to smaller eddies in all directions and the directional nature of the 
primary eddies is gradually lost. Kolmogoroff (1941) concludes that all eddies which are 
much smaller than the primary eddies are statistically independent of them and the 
properties of these small eddies are determined by the local energy dissipation rate per unit 
mass of fluid. For local isotropic turbulence the smallest eddies are responsible for most of 
the energy dissipation and their time scale is given by Kolmogoroff (1941). Turbulence in the 
immediate vicinity of a bubble affects heat and mass transfer rates between the bubble and 
the liquid and may lead to its breakup. 

Kastanek et al. (1993) suggested that the mass transfer in the turbulent bulk-liquid region is 
accomplished by elementary transfer eddies while in the surface layer adjacent to the 
interface turbulence is damped and mass transfer occurs due to molecular diffusion. In 
agreement with the theory of isotropic turbulence, the authors represented the contact time 
as the ratio of the length of elementary transport eddy to its velocity at the boundary 
between the bulk liquid and diffusion layer. Kastanek et al. (1993) argue that the rate of 
mass transfer between the gaseous and the liquid phase is decisively determined by the rate 
of energy dissipation in the liquid phase. 

Kastanek (1977) and Kawase et al. (1987) developed a theoretical model for prediction of 
volumetric mass transfer coefficient in bubble columns. It is based on Higbie's (1935) 
penetration theory and Kolmogoroff's theory of isotropic turbulence. It is believed that 
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turbulence brings up elements of bulk fluid to the free surface where unsteady mass transfer 
occurs for a short time (called exposure or contact time) after which the element returns to 
the bulk and is replaced by another one. The exposure time must either be determined 
experimentally or deduced from physical arguments. Calderbank and Moo-Young (1961) 
and Kawase et al. (1987) developed a correlation relating the rate of energy dissipation to 
turbulent mass transfer coefficient at fixed surfaces. 


1.5 Estimation of gas-phase mass transfer coefficient 

There is a lack of research in the literature on the estimation of the volumetric gas-phase 
mass transfer coefficients kca. On the basis of chemical absorption and vaporization 
experiments Metha and Sharma (1966) correlated the kca values to the molecular diffusivity 
in gas, the superficial gas velocity and static liquid height. Botton et al. (1980) measured kca 
by the chemical absorption method in a SO; (in air)-NazCO3 aqueous solution system in a 
wide range of uc values. Cho and Wakao (1988) carried out experiments on stripping of five 
organic solutes with different Henry's law constants in a batch bubble column with water 
and they proposed two correlations (for single nozzle and for porous tube spargers) for kca 
prediction. Sada et al. (1985) developed also correlation for kca prediction. In the case of 
slurry bubble columns, the authors measured kga by using chemical adsorption of lean CO2 
into NaOH aqueous solutions with suspended Ca(OH)? particles and they developed a 
correlation. Its predictions agree well with those observed by Metha and Sharma (1966) and 
Botton et al. (1980) in a bubble column. 

The gas-phase mass transfer coefficient kc decreases with increasing pressure due to the fact 
that the gas diffusion coefficient is inversely proportional to pressure (Wilkinson et al., 
1992). In the case of bubble columns equipped with single nozzle and porous tube spargers 
the kca value can be calculated by the correlation of Cho and Wakao (1988). 


1.6 Estimation of interfacial area 

The specific gas-liquid interfacial area varies significantly when hydrodynamic conditions 
change. Several methods exist for interfacial area measurements in gas-liquid dispersions. 
These are photographic, light attenuation, ultrasonic attenuation, double-optical probes and 
chemical absorption methods. These methods are effective under certain conditions only. 
For measuring local interfacial areas at high void fractions (more than 20 %) intrusive 
probes (for instance, double optical probe) are indispensable (Kiambi et al., 2001). 
Calderbank (1958) developed a correlation for the specific interfacial area in the case of non- 
spherical bubbles. Akita and Yoshida (1974) derived also their own empirical equation for 
the estimation of the interfacial area. Leonard and Houghton (1963) related the interfacial 
area to the bubble volume using a constant shape factor, which would be 4.84 for spherical 
bubbles. In the case of perforated plates, Kastanek et al. (1993) reported a correlation for the 
estimation of the interfacial area. Very frequently the reliability of estimation of the specific 
interfacial area depends on the accuracy of gas holdup determination. 


2. Effect of bubble shape on mass transfer coefficient 


Deformed bubbles are generally classified as ellipsoids or spherical caps (Griffith, 1960; 
Tadaki and Maeda, 1961). The shapes and paths of larger non-spherical bubbles are 
generally irregular and vary rapidly with time, making the exact theoretical treatment 
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impossible. Bubbles greater than 1.8x10? m in diameter assume mushroom-like or 
spherical-cap shapes and undergo potential flow. Calderbank (1967) argue that the 
eccentricity decreases with increasing viscosity accompanied by the appearance of "tails" 
behind small bubbles and of spherical indentations in the rear surfaces. 

Calderbank et al. (1970) developed a new theory for mass transfer in the bubble wake. Their 
work with aqueous solutions of glycerol covers the bubble size range 0.2-6.0x10 m and 
includes the various bubble shapes as determined by the bubble size and the viscosity of the 
Newtonian liquid. Calderbank and Lochiel (1964) measured the instantaneous mass transfer 
coefficients in the liquid phase for carbon dioxide bubbles rising through a deep pool of 
distilled water. Redfield and Houghton (1965) determined mass transfer coefficients for 
single carbon dioxide bubbles averaged over the whole column using aqueous Newtonian 
solutions of dextrose. Davenport et al. (1967) measured mass transfer coefficients averaged 
over column lengths of up to 3 m for single carbon dioxide bubbles in water aqueous 
solutions of polyvinyl alcohol and ethyl alcohol, respectively. Angelino (1966) has reported 
some shapes and terminal rise velocities for air bubbles in various Newtonian liquids. 
Liquid-phase mass transfer coefficients for small bubbles rising in glycerol have been 
determined by Hammerton and Garner (1954) over bubble diameters ranging from 0.2x102 
m to 0.6x102 m. Barnett et al. (1966) reported the liquid-phase mass transfer coefficients for 
small CO» bubbles (0.5-4.5x105 m) rising through pseudoplastic Newtonian liquids. This 
bubble size range was extended to 3-50x103 m in the data reported by Calderbank (1967). 
Astarita and Apuzzo (1965) presented experimental results on the rising velocity and shapes 
of bubbles in both purely viscous and viscoeleastic non-Newtonian pseudoplastic liquids. 
According to Calderbank et al. (1970) bubble shapes observed in distilled water vary from 
spherical to oblate spheroidal (0.42-1.81x102 m) to spherical cap (1.81-3.79x102 m) with 
increasing bubble size. Over the size range (4.2-70x10? m) the bubbles rise with a zigzag or 
spiral motion and between bubble diameters of 7x10? m (Re=1800) and 18x10? m (Re=5900) 
an irregular ellipsoid shape is adopted and the bubble pulsates about its mean shape. Over 
the bubble size range 1.8-3.0x102 m a transition from irregular ellipsoid to spherical cap 
shape occurs and surface rippling is much more evident. For bubble sizes greater than 3x10- 
? m the bubbles adopt fully developed spherical cap shapes and exhibit little surface 
rippling. These spherical cap bubbles rise rectilinearly. 

Calderbank et al. (1970) developed theory of mass transfer from the rear of spherical-cap 
bubbles. The authors argue that the overall mass transfer coefficients enhance by 
hydrodynamic instabilities in the liquid flow round bubbles near the bubble shape 
transition from spherical cap to oblate spheroid. Calderbank et al. (1970) reported that for 
bubble sizes 1-1.8x102 m a shape transition occurs, the bubble rear surface is gradually 
flattening and becoming slightly concave as the bubble size is increased. The onset of 
skirting is accompanied by a flattening of the bubble rear surface. The authors argue that the 
bubble eccentricity decreases with increasing Newtonian liquid viscosity, though there is a 
tendency towards convergence at large bubble sizes. 

Davidson and Harrison (1963) indicated that the onset of slug flow occurs approximately at 
bubble size/column diameter>0.33. In the case of spherical-cap bubbles it is expected that 
there will be appreciable variations between front and rear surfaces. Behind the spherical- 
cap bubble is formed a torroidal vortex. 

Calderbank et al. (1970) reported that a maximum value of k, occurs shortly before the onset 
of creeping flow conditions and corresponds to a bubble shape transition from spherical cap 
to oblate spheroid. This shape transition and the impending flow regime transition results in 
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instabilities in the liquid flow around the bubble resulting in ki enhancement. The results of 
Zieminski and Raymond (1968) indicate that for CO» bubbles a maximum of kı occurs at 
bubble size of 3x10 m which they attribute to a progressive transition between circulating 
and rigid bubble behavior. 

Calderbank (1967) stated that the theory of mass transfer has to be modified empirically for 
dispersion in a non-isotropic turbulent field where dispersion and coalescence take place in 
different regions. Coalescence is greatly influenced by surfactants, the amount of dispersed 
phase present, the liquid viscosity and the residence time of bubbles. The existing theories 
throw little light on problems of mass transfer in bubble wakes and are only helpful in 
understanding internal circulation within the bubble. The mass-transfer properties of bubble 
swarms in liquids determine the efficiency and dimensions of the bubble column. 

If the viscous or inertial forces do not act equally over the surface of a bubble they may 
cause it to deform and eventually break. A consequence of these dynamic forces acting 
unequally over the surface of the bubble is internal circulation of the fluid within the bubble 
which induces viscous stresses therein. These internal stresses also oppose distortion and 
breakage. 


3. New approach for prediction of gas holdup (Nedeltchev and Schumpe, 
2008) 


Semi-theoretical approaches to quantitatively predict the gas holdup are much more reliable 
and accurate than the approaches based on empirical correlations. In order to estimate the 
mass transfer from bubbles to the surrounding liquid, knowledge of the gas-liquid 
interfacial area is essential. The specific gas-liquid interfacial area, defined as the surface 
area available per unit volume of the dispersion, is related to gas holdup ec and the 
Sauter-mean bubble diameter d; by the following simple relation: 


z Beg 1 
7 (1) 
Strictly speaking, Eq. (1) (especially the numerical coefficient 6) is valid only for spherical 
bubbles (Schügerl et al., 1977). 
The formula for calculation of the interfacial area depends on the bubble shape. Excellent 
diagrams for bubble shape determination are available in the books of Clift et al. (1978) and 
Fan and Tsuchiya (1990) in the form of log-log plots of the bubble Reynolds number Res vs. 
the Eótvós number Eo with due consideration of the Morton number Mo. A comparison 
among the experimental conditions used in our work and the above-mentioned standard 
plots reveals that the formed bubbles are no longer spherical but oblate ellipsoidal and 
follow a zigzag upward path as they rise. Vortex formation in the wake of the bubbles is 
also observed. The specific interfacial area a of such ellipsoidal bubbles is a function of the 
number of bubbles Ns, the bubble surface Sg and the total dispersion volume Vootal 
(Painmanakul et al., 2005; Nedeltchev et al., 2006a, b, 20072): 


a= N55 = Nps (2) 
Viotal AH 


where A denotes the column cross-sectional area. The number of bubbles Ng can be deduced 
from the bubble formation frequency fg and bubble residence time (Painmanakul et al., 2005): 
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Ho Uan (3) 


N.= 
B f TA Up 


where Qc is the volumetric gas flow rate, ug is the bubble rise velocity and Vg is the bubble 
volume. The substitution of Eq. (3) into Eq. (2) yields: 


_ fsg 
a= [os (4) 


The bubble rise velocity ug can be estimated from Mendelson's (1967) correlation: 


t= | T + (5) 
Pide 2 


This equation is particularly suitable for the case of ellipsoidal bubbles. 
The volume of spherical or ellipsoidal bubbles can be estimated as follows: 


3 2 
Wwe h (6) 


6 3 


If some dimensionless correction factor f. due to the bubble shape differences is introduced, 
then Eqs. (1) and (4) might be considered equivalent: 


66, fas 


= 7 
4. fe ros (7) 
Rearrangement of Eq. (7) yields: 
d, fS 
= SJ B~B 8 
EG Js 6Au, ( ) 


The surface Ss of an ellipsoidal bubble can be calculated as follows (Nedeltchev et al., 2006a, 


b, 20072): 
2 2 
1 
m Lo Vik +?) (9) 
2 I) 2e (1-8) 
where e is the bubble eccentricity. It can be calculated as follows: 


" - (t (92) 


An oblate ellipsoidal bubble is characterized by its length 1 (major axis of the ellipsoid) and 
its height h (minor axis of the ellipsoid). The ellipsoidal bubble length 1 and height h can be 
estimated by the formulas derived by Tadaki and Maeda (1961) and Terasaka et al. (2004): 
for 2«Ta«6: 


d 
s - 10a 
Line (10a) 
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h -1.3d, Ta 0??? (10b) 
for 6<Ta<16.5: 
d 
BET (a) 
h -1.85d,Ta ^6 (11b) 
where 
Ta = Reg Mo??? (12) 
Reg = s (13) 
L 
ut 
Moa E (14) 
PLOL 


It is worth noting that the major axis of a rising oblate ellipsoidal bubble is not always 
horizontally oriented (Yamashita et al., 1979). The same holds for the minor axis of a rising 
oblate ellipsoidal bubble, i.e. it is not necessarily vertically oriented (Akita and Yoshida, 1974). 
Equations (10a)-(14) were used to calculate both | and h values under the operating 
conditions examined. The Morton number Mo is the ratio of viscosity force to the surface 
tension force. The Tadaki number Ta characterizes the extent of bubble deformation; the Ta 
values fell always in one of the ranges specified above. This fact can be regarded as an 
additional evidence that the bubbles formed under the operating conditions examined are 
really ellipsoidal. 

The above correlations (Equations (10a)—(14)) imply that one needs to know a priori the 
bubble equivalent diameter de. Very often in the literature is assumed that de can be 
approximated by the Sauter-mean bubble diameter d.. The latter was estimated by means of 
the correlation of Wilkinson et al. (1994): 


E —0.04 3 -0.12 0.22 
gpa, L. 8.8 Uc A. A (15) 
OL OL 8H, Po 


Equation (15) implies that the bubble size decreases as the superficial gas velocity uc or the 
gas density pc (operating pressure P) increase. The calculated d, values for all liquids 
examined imply an ellipsoidal shape. Equation (5) along with Eq. (15) (for ds estimation) 
was used also to calculate the bubble Reynolds number Res (Eq. (13)) needed for the 
estimation of both | and h values. 

The bubble equivalent diameter d. of an ellipsoidal bubble can be also calculated from Eq. 
(6) by assuming a sphere of equal volume to the volume of the ellipsoidal bubble: 


d, - (Ph 


) 1/3 (16) 
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Estimating the characteristic length of ellipsoidal bubbles with the same surface-to-volume 
ratio (the same d, value as calculated from Eq. (15)) required an iterative procedure but led 
to only insignificantly different values than simply identifying the equivalent diameter de 
with ds when applying Eqs. (10a—b) or (11a-b). In other words, the differences between 
bubble diameters estimated by Eq. (15) and Eq. (16) are negligibly small. 


Liquid D. Gas Gases P PL pL OL 
4 [m] Sparger Used [MPa] [kg-m?] [103 Pa.s] [103 N-m'1] 

Acetone 0.095 D4 air 0.1 790 0.327 23.1 
Anilin 0.095 D4 air 0.1 1022 4.4 43.5 
Benzene 0.095 D4 air 0.1 879 0.653 28.7 

0.0055 D1, D2 ; 
1-Butanol 0.102 DA Nz air, He 0.1-4.0 809 2.94 24.6 

. air, He, Ho, 
Carbon tetrachloride ^ 0.095 D4 CO 0.1 1593 0.984 26.1 
2 

Cyclohexane 0.095 D4 air 0.1 778 0.977 24.8 
Decalin 0102 D1, D2 N, He 0.1-4.0 884 2.66 32.5 
1,2-Dichloroethane 0.095 D4 air 0.1 1234 0.82 29.7 
1,4-Dioxane 0.095 D4 air 0.1 1033 1.303 32.2 
Ethanol (96 %) 0.102 D1, D2, D3 N,He 0.1-4.0 793 1.24 22.1 
Ethanol (99 %) 0.095 D4 air 0.1 791 1.19 22.1 
Ethyl acetate 0.095 D4 air 0.1 900 0.461 23.5 
Ethyl benzene 0.095 D4 air 0.1 867 0.669 28.6 
Ethylene glycol os D1, D2, DA Nz, air, He 0.1-4.0 1112 19.9 47.7 
Gasoline 0.102 D1 N2 0.1-4.0 692 0.464 21.6 
LigroinA (b. p. 90-110 °C) 0.095 D4 air 0.1 714 0.470 20.4 
Ligroin B(b. p. 100-140°C) 0.095 D4 air 0.1 729 0.538 21.4 
Methanol 0.095 D4 air 0.1 790 0.586 22.2 
Nitrobenzene 0.095 D4 air 0.1 1203 2.02 38.1 
2-Propanol 0.095 D4 air 0.1 785 2.42 21.1 

0.095 D1, D2 . 
Tap water 0102 D3, D4 Nz air, He 0.1-4.0 1000 1.01 72.74 
Tetralin 0.095 D4 No ai, 01 968 2418 34.9 

CO» 

0.095 ; 
Toluene 0.102 D1, D2, D4 Nz air, He 0.1-4.0 866 0.58 28.5 
Xylene 0.095 pa Nean He 91 863 0.63 28.4 


Hə, CO; 


Table 1. Properties of the organic liquids and tap water (293.2 K) 
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Our semi-theoretical approach is focused on the derivation of a correlation for the 
correction term f. introduced in Eq. (7). Many liquids covering a large spectrum of 
physicochemical properties, different gas distributor layouts and different gases at 
operating pressures up to 4 MPa are considered (Nedeltchev et al., 2007a). As many as 386 
experimental gas holdups were obtained in two bubble columns. The first stainless steel 
column (D.-0.102 m, Ho-1.3 m) was equipped with three different gas distributors: 
perforated plate, 19 x Ø 1x10? m (D1), single hole, 1 x Ø 4.3x10? m (D2) and single hole, 
1 x Ø 1x10? m (D3) (Jordan and Schumpe, 2001). In the second plexiglass column 
(D-=0.095 m, Ho-0.85 m) the gas was always introduced through a single tube of 3x10? m in 
ID (D4) (Öztürk et al., 1987). The ec values were measured in 21 organic liquids, 17 liquid 
mixtures and tap water (see Tables 1 and 2). 


Key D. Gas Gas P PL PL OL 


Liquid Mixture Fig.9 [m] Sparger Used [MPa] [kg:m3] [10-3 Pa-s] [103 N-m-1] 


Benzene/ Cyclohexane 6.7 % 0.095 D4 air 01 865 0.634 27.6 
Benzene/ Cyclohexane 13.4 % 0.095 D4 air 01 854 0.628 26.9 
Benzene/ Cyclohexane 31.5 % 0.095 D4 air 01 834 0.631 26.2 
Benzene/ Cyclohexane 54 % 0.095 D4 air 01 814 0.672 254 
Benzene/ Cyclohexane 78.5 % 0.095 D4 air 01 797 0.772 24.9 
Benzene/ Cyclohexane 90 % 0.0050 D4 air 01 787 0.858 24.9 
Glycol 22.4 % / Water 0.05 D4 air 01 1043 2.32 53.8 
Glycol 60 % /Water 0.095 D4 air 0.1 1072 5.6 52.0 
Glycol 80 % / Water 0.05 D4 air 01 1091 9.65 51.0 


Toluene/ Ethanol 5 % 
Toluene/ Ethanol 13.6 % 
Toluene/ Ethanol 28 % 
Toluene/ Ethanol 55 % 
Toluene/ Ethanol 73.5 % 
Toluene/ Ethanol 88.5 % 
Toluene/ Ethanol 94.3 % 
Toluene/ Ethanol 97.2 % 


0.000 D4 air 04 863 0.578 27.6 
0.095 D4 air 04 859 0.587 27.3 
0.000 D4 air 01 852 0.616 25.5 
0.095 D4 air 01 838 0.731 25.0 
0.000 D4 air 01 823 0.961 24.2 
0.000 D4 air 01 807 1.013 23.3 
0.095 D4 air 01 800 1.103 22.7 
0.095 D4 air 01 796 1.135 22.2 


O.-o0oZzznv nm—^-—-iugüomnumnocnus»s 


Table 2. Properties of the liquid mixtures (293.2 K) 


In both tables are listed the different combinations of liquids, gases, gas distributors and 
operating pressures that have been used. It is worth noting that in a 0.095 m in ID bubble 
column equipped with a sparger D4 every liquid or liquid mixture was aerated with air. 
Table 1 shows that in the case of few liquids (carbon tetrachloride, tetralin, toluene and 
xylene) some other gases have been used. It should be mentioned that in the case of 0.102 m 
in ID bubble column no air was used (only nitrogen and helium). 

The gas holdups ec in 1-butanol, ethanol (96 76), decalin, toluene, gasoline, ethylene glycol and 
tap water were recorded by means of differential pressure transducers in the 0.102 m stainless 
steel bubble column operated at pressures up to 4 MPa. The following relationship was used: 
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"E AP 6 gas AP sas (17) 
AB, gas 


where AP is the pressure difference between the readings of both lower (at 0 m) and upper 
(at 1.2 m) pressure transducers. The subscript “no gas" denotes the pressure difference at 
the clear liquid height Ho, whereas the subscript “gas” denotes the pressure difference at the 
aerated liquid height H. The gas holdups sc in all other liquids and liquid mixtures were 
estimated by visually observing the dispersion height under ambient pressure in the 0.095 m 
in ID bubble column. The upper limit (transitional gas velocity) of the homogeneous regime 
(transition gas velocity Utrans) was estimated by the formulas of Reilly et al. (1994). 

Most of the gas-liquid systems given in Tables 1 and 2 were characterized with Tadaki 
numbers Ta lower than 6 and thus Eqs. (10a—b) for the estimation of both bubble length 1 
and bubble height h were applied. Only in the case of ethylene glycol-(helium, air and 
nitrogen), 1-butanol-(helium and air) and decalin-helium the Ta values exceeded 6 and 
then Eqs. (11a-b) were used. 

It was found that the dimensionless correction factor f. can be correlated successfully to both 
the Eótvós number Eo and a dimensionless gas density ratio: 


0.07 ( _ ) d? -0.22 0.07 
f, -078Eo 2 a -078| 85-524 es =) (2) (18) 
Pa OL 1.2 


G 


where pc'*f is the reference gas density (1.2 kg-m? for air at ambient conditions: 293.2 K and 
0.1 MPa). All experimental gas holdup data (386 points) were fitted with an average error of 
9.675. The dimensionless gas density ratio is probably needed because the correlation of 
Wilkinson et al. (1994) was derived for pressures up to 1.5 MPa only, whereas the present 
data extend up to the pressure of 4 MPa. It is worth noting that Krishna (2000) also used 
such a dimensionless gas density ratio for correcting his correlations for large bubble rise 
velocity and dense-phase gas holdup. 

Figure 1 illustrates the decrease of the product f«(pc/1.2)997 with increasing Eo. At smaller 
bubble sizes (with shapes approaching spheres) Eo will be lower and thus f. higher 
(gradually approaching unity). It is worth noting that most of the liquids are characterized 
with Eo values in a narrow range between 2 and 8. 

Figure 2 illustrates that the correction factor f. increases with gas density pc (operating 
pressure) leading to bubble shrinkage. For example, the correction term f. decreases in the 
following sequence: toluene > ethanol > decalin > 1-butanol > ethylene glycol. The smallest 
bubble size is formed in toluene, whereas the largest bubble size is formed in the case of 
ethylene glycol. When very small (spherical) bubbles are formed, the correction factor f. 
should be equal to unity and both expressions for the interfacial area should become 
identical (see Eq. (7)). 

Figures 3 and 4 exhibit that the experimental gas holdups £c measured in 1-butanol, decalin 
and toluene at pressures up to 4 MPa can be predicted reasonably well irrespective of the 
gas distributor type. 

The same result holds for ethylene glycol and tap water (see Fig. 5). The successful prediction 
of gas holdups in ethylene glycol should be regarded as one of the most important merits of 
the presented method since the viscosity is much higher than that of the other liquids. 
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Fig. 1. Product f(pc/1.2)907 as a function of Eo for 12 organic liquids, 2 liquid mixtures and 
tap water at ambient pressure (Gas distributor: D4 unless specified in the legend) 
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Fig. 2. Correction factor f. as a function of gas density pc in five organic liquids (Gas 
distributor: D1; D.: 0.102 m) 
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Fig. 3. Parity plot for gas holdups in 1-butanol and decalin sparged with nitrogen through 
gas distributors D1 and D2 at various pressures 
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Fig. 4. Parity plot for gas holdups in toluene and decalin sparged with nitrogen through gas 
distributors D1 and D2 at various pressures 
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Fig. 5. Parity plot for gas holdups in ethylene glycol and tap water aerated with nitrogen 
through gas distributors D1 and D2 at various pressures 
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Fig. 6. Parity plot for gas holdups in ethanol (96 76) sparged with nitrogen through gas 
distributors D1—D3 at various pressures 
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Figure 6 shows, for ethanol (96 76) as an example, that the gas distributor type is not so 
important. The same holds for 1-butanol, decalin and toluene (Figs. 3-5) and tap water. This 
fact is in agreement with the work of Wilkinson et al. (1992) who stated that once the hole 
size of the gas distributor is greater than 1-2x10? m, then it has no significant effect on the 
gas holdup. 

Eight organic liquids and tap water have been aerated not only with air or nitrogen but also 
with other gases (helium, hydrogen and carbon dioxide). Figure 7 shows that the developed 
model is capable of predicting satisfactorily the experimental gas holdups at these operating 
conditions. 


@ Carbon tetrachloride/CO2 
© Carbon tetrachloride/He 
A Carbon tetrachloride/H2 
94 Tetralin/CO2 
$$ Xylene/CO2 
El Xylene/He 
O Xylene/H2 
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A 1-Butanol/He (D2) 
+ Ethanol (96 %)/He (D2) 
x Ethanol (96 %)/He (D3) 
O Toluene/He (D1) 
@ Toluene/He (D2) 
@ Decalin/He (D1) 
A Ethylene glycol/He (D2) 
E Water/He (D1) 
0.01 0.03 0.05 0.07 0.09 0.11 0.13 X Water/He (D2) 

£c (calculated) [-] 2? Water/He (D3) 


0.07 


0.05 


£g; (experimental) [-] 


0.03 


Fig. 7. Parity plot for gas holdups ec in 8 organic liquids and tap water aerated with other 
gases (helium, hydrogen and/or carbon dioxide) at ambient pressure. Gas distributor: D4 
unless specified in the legend 


Figure 8 shows that the model predicts reasonably well the experimental gas holdups 
measured in 15 organic liquids at ambient pressure. This fact should be regarded as further 
evidence that by the introduction of a correction term the presented method becomes 
generally applicable. 

Table 2 and Figs. 3-8 reveal that our approach is applicable not only to tap water but also to 
organic liquids covering the following ranges of the main physicochemical properties: 
692x(pi / kg 13)x1593, 0.327x10?x(u;/Pa s)x 19.9x105, 20.4x103x(o1/ N m1)x47.7x103. 
Figure 9 exhibits that the proposed method for gas holdup prediction along with the new 
correction factor (Eq. (18)) is also applicable to various liquid mixtures. Table 2 shows that 
the examined liquid mixtures cover the following ranges of the main physicochemical 
properties: 787<(pL/ kg m)x1091, 0.578x103x(qu / Pa s)x9.65x103 and 22.2x103x(o1/ N : md) 
«53.8x103. 
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Fig. 8. Parity plot for gas holdups in 15 organic liquids aerated with air by means of gas 
distributor D4 (D.=0.095 m) at ambient pressure 


0.13 @ Mixture A 
Mixture B 
A Mixture C 
Mixture D 
a= Mixture E 
M Mixture F 

Mixture G 
Mixture H 
Mixture I 

Mixture J 

Mixture K 
Mixture L 
Mixture M 
Mixture N 
Mixture O 
Mixture P 

Mixture Q 


0.11 


o 
> 
S 
+P 0 


£g; (experimental) [-] 
x] 


xE >xo0 


0.01 0.03 0.05 0.07 0.09 0.11 0.13 
£c (calculated) [-] 


Fig. 9. Parity plot for gas holdups in 17 liquid mixtures sparged with air (gas distributor D4, 
D.-0.095 m) at ambient pressure. The legend keys are given in Table 2 
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Nedeltchev and Schumpe (2008) have shown that in the homogeneous flow regime the 
semi-theoretical approach improves the gas holdup predictions and turns out to be the 
most reliable one (since it produces better gas holdup predictions than the empirical 
correlations of both Hammer et al. (1984) and Akita and Yoshida (1974)). It is worth 
underlying that our approach is applicable at ucSugas, where the transition gas velocity 
Utrans can be predicted by the correlations of Reilly et al. (1994). In this work, the Utrans values 
were always less than 0.04 m:s-. 


4. New approach for prediction of volumetric liquid-phase mass transfer 
coefficient (Nedeltchev et al., 2006a, b, 2007a) 


The Kia value can be predicted if one knows how to estimate both the liquid-side mass 
transfer coefficient kų and the interfacial area a. These parameters depend on the bubble 
diameter. Both kų and a are closely related to the bubble bed hydrodynamics. k 
incorporates the effects of the complex liquid flow field surrounding the rising gas bubbles. 
The interfacial area a inherently reflects the system bubble behavior. 

Gas bubbles are nonspherical, except when interfacial tension and viscous forces are much 
more important than inertial forces (Clift et al., 1978). Most of the bubble sizes in this work fall 
in the range of 1.4-6x103 m. Such bubbles are no longer spherical and follow a zigzag or 
helical upward path. Viscous drag is augmented by vortex formation in the wake, and velocity 
of rise remains fairly constant over the bubble size range (Miller, 1974). 

As a rough approximation, Higbie (1935) assumed that the average time of exposure for mass 
transfer (called contact time or exposure time for mass transfer) can be estimated as follows: 


Bubble surface 


c : (19) 
Rate of surface formation 
In our approach all theoretical kia calculations were based on this definition of the 
gas-liquid contact time tc. The latter characterizes the residence time of the micro eddies 
(responsible for mass transfer in the liquid film) at the interface, i.e. at the bubble surface. It 
is practically impossible to measure directly the contact time t, so usually it is estimated 
from correlations for the bubble surface and the rate of surface formation. The contact time 
is the time required for the bubble to rise through its equivalent spherical diameter. Higbie's 
(1935) model predicts a decrease in kų with increase in diameter when in fact it increases 
quite rapidly. Higbie's theory does not predict an effect of distance from the point of release. 
The simple model of Higbie (1935) in which the contact time is approximated by macro- 
scale parameters fails for turbulent bubbling conditions where bubble swarms exist. 
Higbie's (1935) model assumes direct contact between the gas phase and the bulk liquid. 

The contact time takes part in the evaluation of the liquid-phase mass transfer coefficient: 


Ee (20) 


The equation is valid explicitly for rigid spherical bubbles. Higbie (1935) postulated that a 
gas bubble moving through a liquid splits the liquid at its advancing tip. The penetration 
theory assumes unsteady-state absorption of a gas by a fluid element adjacent to the 
surface. The element moves at a uniform velocity from the front of the bubble to the rear as 
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penetration into it occurs. Timson and Dunn (1960) showed that the bubble surface increases 
when the spherical bubble is deformed into an ellipsoid. This leads to higher contact time 
and thus lower ki value. 

In order to calculate the volumetric liquid-side mass transfer coefficient kra one also needs 
to know how to calculate the interfacial area a. The formula for its calculation depends on 
the bubble shape. The specific interfacial area a is a function of the number of bubbles Ns, 
the bubble surface area Sg and the total dispersion volume Vio (see Eq. (2)). 

The theoretically calculated kia values were obtained as a product of both Eq. (20) (using 
Eq. (19) for the contact time) and Eq. (2). Experimentally, on the other hand, it is much easier 
to measure the product of kų and a than the individual values. In the homogeneous flow 
regime there is a narrow bubble size distribution and the researchers use frequently the 
Sauter-mean bubble diameter d, for their mass transfer calculations. The bubbles rise 
uniformly in nearly straight lines and have roughly uniform bubble size (Krishna, 2000; 
Lucas et al., 2005). Therefore, it is a reasonable simplification in the homogeneous flow 
regime to use the mean bubble diameter and disregard the bubble size distribution. 

For the sake of theoretical prediction of the kra values by means of Eqs. (2), (19) and (20) one 
needs to calculate the bubble surface Sz, the rate of surface formation Rs and the number of 
bubbles Ns. The calculation of the bubble surface Sg depends on the bubble shape 
(Painmanakul et al., 2005). Two excellent diagrams for bubble shape determination are 
available in the books of Clift et al. (1978) and Fan and Tsuchiya (1990) in the form of 
log-log plots of the bubble Reynolds number Res vs. the Eótvós number Eo with due 
consideration of the Morton number Mo. A comparison among the experimental conditions 
used in this work and the above-mentioned standard plots reveals that the formed bubbles 
are no longer spherical but oblate ellipsoidals that follow a zigzag upward path as they rise. 
Vortex formation in the wake of the bubbles is also observed. 

An oblate ellipsoidal bubble is characterized by its length 1 (major axis of the ellipsoid) and 
its height h (minor axis of the ellipsoid). The surface area Sg of an ellipsoidal bubble is 
calculated on the basis of Eqs. (9) and (9a) (Fan and Tsuchiya 1990; Nedeltchev et al., 2006a, 
b). In order to calculate both ellipsoidal bubble length 1 and height h, the correlations of 
Terasaka et al. (2004) (for 2«Ta«6) were used (see Eqs. (10a-b)). The Tadaki numbers Ta fell 
always in the range specified above. However, the correlations of Terasaka et al. (2004) 
imply that one needs to know a priori the bubble equivalent diameter de. 

Very often in the literature is assumed that de can be approximated by the Sauter-mean 
bubble diameter ds. Bubble shape, motion and any tendency for the interface to ripple, 
fluctuate or otherwise deform are all related to bubble size. In turn, bubble size is 
determined by the physical characteristics of the system and the operating conditions. The 
d, value was estimated by means of the correlation of Wilkinson et al. (1994) which is one of 
the most frequently cited in the literature. This equation implies that the bubble size 
decreases as superficial gas velocity uc or gas density pc (operating pressure P) increase. 
The calculated d; values for all liquids examined correspond to an ellipsoidal shape. 

The bubble equivalent diameter d. of an ellipsoidal bubble was calculated by assuming a 
sphere of volume equal to the volume of the ellipsoidal bubble (see Eqs. (6) and (16)). 
Estimating the characteristic length of ellipsoidal bubbles with the same surface-to-volume 
ratio (the same d; value as calculated from Wilkinson et al.'s (1994) correlation) required an 
iterative procedure but led to only insignificantly different values than simply identifying 
the equivalent diameter de with ds when applying the correlations for both | and h (see 
Terasaka et al., 2004). 
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The equivalent bubble diameter de is needed also for the calculation of the bubble rise velocity 
up. The latter was estimated by means of the correlation of Mendelson (1967). This equation is 
particularly suitable for the case of ellipsoidal bubbles. Mendelson's (1967) correlation for ug 
along with Wilkinson et al.'s (1994) relationship for d; estimation were also used to calculate 
the bubble Reynolds number Res (Eq. (13)) needed for the estimation of both 1 and h values. 
The bubble rise velocity ug and both the bubble length 1 and height h of an ellipsoidal 
bubble take part in the calculation of the rate of surface formation Rs (Higbie, 1935). For 
oblate ellipsoidal bubbles (see Nedeltchev et al., 2006a, b): 


Pen (I-hy 
n EU e, 2 


Equation (21) is needed for the calculation of the contact time t. (Eq. 19) and thus the 
liquid—side mass transfer coefficient ki, (Eq. 20). 

The number of bubbles Ns was deduced from the bubble residence time (aerated liquid 
height H divided by the bubble rise velocity ug) and the bubble formation frequency fg 
(Painmanakul et al, 2005; Nedeltchev et al, 2006a, b). The latter was expressed as 
volumetric flow rate Oc divided by the bubble volume Vg (see Eqs. (3) and (6)). The 
interfacial area can be calculated as a function of bubble frequency fs, bubble surface Ss, 
bubble rise velocity ug and cross section (see Eq. (4)). Therefore, the theoretical kya values 
for ellipsoidal bubbles can be calculated as a product of ki, value (estimated by means of 
Eqs. (5), (9), (10a-b), (20) and (21)) and a value (estimated by means of Eq. (4)). 

The theoretical kra values have to be multiplied with some correction factor f. in order to 
match satisfactorily the experimental kra values. Specifically, the accurate kra values for 
ellipsoidal bubbles should be calculated by means of the following formula: 


4D, Ry feSp 
mSp Aug 


ka= f, (22) 


The introduction of the correction factor f. could be attributed to the fact that even when Eq. 
(19) is modified for oblate ellipsoidal bubbles it does not take into account the effect of the 
bubble wake and surface disturbances. Due to these supplementary effects on the mass 
transfer rate some additional term should be introduced. 

Nedeltchev et al. (2006a) have illustrated that at pressures up to 1 MPa a good prediction of 


the experimental kj a values is achieved when the correction factor is expressed as a function 
only of the Eótvós number Eo: 


f, = 0.185E0°”*” (23) 


where 


gle- Pc) a 
OT 


Eo= (24) 


The Eötvös number Eo represents the gravitational force-to-surface tension force ratio. As 
mentioned above, the bubble shape depends on the Eo value. Eqs. (23) and (24) show that as 
pc increases and especially as bubble size reduces, Eo and thus f. become lower. 
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The above-described theoretical method was applied to conduct a systematic comparison 
between predicted and experimental kya results over a wide range of physicochemical 
properties, which are characteristic for gas-liquid systems. Both pure organic liquids and 
liquid mixtures as well as tap water were considered. In such a way, the capability of the 
correction method to predict the numerous experimental kra data available in the literature 
was assessed. Following the introduction of an additional correction term accounting for the 
gas density effect, the predicted kra values at high pressures were improved by as many as 
18 % (in terms of the maximum relative error). 

Volumetric liquid-phase mass transfer coefficients kta measured in 18 organic liquids, 14 
liquid mixtures and tap water at 293. 2 K (Tables 1 and 2) have been considered. Results for 
1-butanol, ethanol (96 476), toluene, decalin and tap water were reported by Jordan and 
Schumpe (2001) for pressures P of 0.1, 0.2, 0.5, 1, 2 and 4 MPa, respectively. The kra data in 
the other pure and organic liquids were measured at ambient pressure by Oztürk et al. 
(1987). These authors presented also some kra data in tap water. All kra data in the present 
approach refer to the dispersion volume. 

The bubble column used by Jordan and Schumpe (2001) had an inner diameter of 0.1 m and 
a height of 2.4 m. Three different gas distributors were used: perforated plate, 19 x Ø 1x10 
m (D1), single hole, 1 x Ø 4.3x10? m (D2) and single hole, 1 x Ø 1x10? m (D3). The clear 
liquid height was set at 1.3 m. Jordan and Schumpe (2001) measured their kra values with a 
dynamic oxygen desorption technique. The oxygen desorption from the liquid was traced 
with an optical probe (MOPS, Comte, GmbH, Hannover, Germany) based on fluorescence 
extinction by oxygen. By dissolving the fluorophore in the liquid, rather than fixing it to the 
tip of the glass fiber, an instantaneous sensor response was achieved. It is worth noting that 
the Dr, values for oxygen were either measured (Öztürk et al., 1987) or calculated by the 
correlation of Schumpe and Lühring (1990). 

Öztürk et al. (1987) carried out their kra experiments in a jacketed glass bubble column of 
0.095 m in ID. The clear liquid height was set at 0.85 m. A single tube of O 3x10? m ID was 
used as the gas distributor (D4). Air, nitrogen, hydrogen or helium were employed as the 
gas phase. 

The kra values were measured by dynamic oxygen absorption or desorption methods. The 
oxygen fugacity in the liquids was measured with a polarographic oxygen electrode 
(WTW-EO 90) inserted horizontally at half the dispersion height. The electrode response 
time was 3 s in water and less in most organic liquids. For absorption runs, oxygen was 
desorbed by sparging nitrogen. After disengagement of all nitrogen bubbles, a preadjusted 
air flow was fed by switching two magnet valves, and the increase in oxygen fugacity was 
recorded. For desorption runs, oxygen-free inert gas (nitrogen, hydrogen or helium) was 
sparged into air-saturated liquid. 

In the present work only kia data obtained in the homogeneous flow regime of bubble 
column operation were analyzed. The upper boundary (so-called transitional gas velocity) 
of this flow regime was determined by the formulas of Reilly et al. (1994). According to the 
bubble shape diagrams presented by Clift et al. (1978) and Fan and Tsuchiya (1990) the 
bubbles formed under all operating conditions were oblate ellipsoidals. 

By means of a nonlinear regressional routine applied to the experimental k;a data the 
following new correction factor was derived: 
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0.15 (p = yd 0.94 0.15 
f, =0.124E0° | PG. | =0.124| EPL Pete | | PG (25) 
c ref 


where pc'*f is the reference gas density (1.2 kg m? for air at ambient conditions: 293.2 K and 
0.1 MPa). 263 experimental kra values are fitted with an average relative error of 10.4 %. It is 
worth noting that this correction factor fe combines the individual corrections of both kų and 
a due to the ellipsoidal shape of bubbles. The average relative error without the gas density 
correction term is 14.9 %. 

The dimensionless gas density ratio is the additional correction term that has not been 
reported hitherto. Its introduction is needed because the correlation of Wilkinson et al. 
(1994) was tested only up to 1.5 MPa. In our work we correlate k; a data at pressures of 2 and 
4 MPa that could not be fitted without this term. Moreover, the correlations of Terasaka et 
al. (2004) have been derived under ambient pressure. To the best of our knowledge, these 
equations have not been validated under high pressure. It is worth noting that such a 
dimensionless gas density correction factor has been also used by Krishna (2000) for correcting 
his correlations for large bubble rise velocity and dense-phase gas holdup. The different 
correction factors reported in our previous papers (Nedeltchev et al., 2006a, b) are due to both 
lower number of liquids considered and lower operating pressures (up to 1 MPa). 

Figs. 10a-10c show the parity plots of kta values measured in four organic liquids aerated 
with three different distributors (D1, D2 and D3) at pressures up to 4 MPa. 

It is clear that by means of the new correction factor (Eq. (25)) the kra values can be 
predicted reasonably well (within + 20 76). It is worth noting that the kia values obtained in 
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Fig. 10a. Parity plot of kra values measured in 1-butanol and decalin at pressures up to 4 MPa 
(Gas distributors: D1 and D2) 
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Fig. 10b. Parity plot of kra values measured in toluene and decalin at pressures up to 4 MPa 
(Gas distributors: D1 and D2) 
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Fig. 10c. Parity plot of kra values measured in ethanol (96 %) at pressures up to 4 MPa 
(Gas distributors: D1, D2 and D3) 
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tap water aerated with gas spargers D1 and D4 can also be predicted satisfactorily by means 
of Eq. (25). One of the merits of our work is the successful prediction of kia data at very high 
pressure (up to 4 MPa). The presented theoretical approach predicts satisfactorily the kga 
values measured by making use of four different gas distributors (D1-D4). It is known that 
the sparger has no substantial influence on both the gas holdup and mass transfer, provided 
that the diameter of the sparger holes is not too small (Wilkinson et al., 1994). 

Fig. 11a exhibits that the correction method predicts satisfactorily the kta values measured 
in 15 different organic liquids aerated with gas sparger D4 at ambient pressure. In the case 
of xylene and tetralin not only air (nitrogen) but also helium (He) and hydrogen were used. 
Ethanol (96 %) was also sparged with helium. 

The new correction method applies not only to pure organic liquids but also to liquid 
mixtures. Fig. 11b shows that the predicted kra values are in reasonable agreement with the 
experimental results. The keys for the different organic mixtures used are given in Table 2. 
All 14 organic mixtures were aerated with a gas sparger D4. As in some other parity plots, 
high kra values tend to be overpredicted. This indicates the onset of coalescence at high gas 
velocity (flow regime transition). In Fig. 2b are included also predicted and experimental kia 
values obtained in ethanol (96 %), toluene and 1-butanol aerated with helium (He) at 
different pressures by using different gas spargers (D2 and D3). These data demonstrate 
again that the developed theoretical approach along with the new correction factor (Eq. (25)) 
is valuable and capable of predicting numerous experimental kra data. 
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Fig. 11a. Parity plot of kra values measured in 15 different organic liquids at ambient 
pressure (Gas distributor: D4) 
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Fig. 11b. Parity plot of kra values measured in 14 different liquid mixtures (P=0.1 MPa, 
sparger D4) as well as in ethanol, toluene and 1-butanol aerated with helium (He) at 
different pressures and gas spargers (D2, D3). The legend keys for the liquid mixtures are 
given in Table 2. 
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Fig. 12 shows the product f.(pc/1.2)945 (see Eq. (25)) as a function of Eo. This product 
increases with Eo for 13 organic liquids at ambient pressure. The trend holds also for 
different gases (air, nitrogen, helium and hydrogen). The f.(pc/1.2)?15 values fall in the 
range of 0.3-1.1 which corresponds to Eo values from 2 to 10. It is worth noting that the f. 
values for the liquids and operating conditions examined in this work are always lower than 
unity (0.3—0.8). In the case of the liquid mixtures specified in Table 2 the f.(pc/1.2)915 vs. Eo 
relationship follows the same trendline as the one shown in Fig. 12. For the sake of brevity 
this figure is omitted. 

Nedeltchev et al. (2006a) have shown that the values of the correction factor (Eq. (23)) fall in 
line with the predictions of the correlation of Wellek et al. (1966) for the bubble aspect ratio E: 


E—1 - 0.163Eo9977 (26) 

where 
E=— 27 
? (27) 


This similarity holds also for Eq. (25). Figure 13 exhibits that the f.(p¢/1.2)°5 trendline at 
Eox5 is very close to the predictions of Eq. (26). Beyond this Eo value the values of the 
product f.(pc/1.2)945 become systematically higher than (E-1) values. 
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Fig. 12. Product f.(pc¢/1.2)°5 as a function of the Eötvös number Eo for 13 different organic 
liquids at ambient pressure (Gas distributor: D4 unless otherwise mentioned) 


It was found that the correction factors f. gradually decrease with the increase of gas density 
pc (operating pressure). As the superficial gas velocity uc increases, a gradual f. reduction is 
also observed. Both trends are illustrated for four organic liquids in Fig. 13. The decreasing 
f. trends are explainable in terms of Eq. (25). The new correction factor f. is primarily 
dependent on the Eótvós number Eo which involves the bubble equivalent diameter de. As 
Uc or pc increases, the bubble size decreases (according to the correlation of Wilkinson et al. 
(1994)) leading to lower Eo and thus lower f. (this trend is stronger than the effect of the gas 
density ratio). 

Larger ellipsoidal bubbles are characterized with higher f. values since larger wakes or 
vortices are formed behind them which enhance the mass transfer coefficient (see Lochiel 
and Calderbank, 1964; Miller, 1974) and thus the correction factor f. becomes closer to unity. 
In addition, with larger ellipsoidal bubbles more surface disturbances occur which also 
enhance the mass transfer characteristics. Fan and Tsuchiya (1990) argue that larger 
ellipsoidal bubbles (larger Res) will have larger shedding frequency of vortex pairs from a 
bubble which will increase the mass transfer coefficient across the base of an ellipsoidal 
bubble and thus the overall liquid—phase mass transfer coefficient. It is worth noting that 
Eq. (20) is valid at very high bubble Reynolds numbers Res (Calderbank, 1967). It means that 
in the same liquid larger ellipsoidal bubbles will have higher Res and thus their correction 
factors f. will be higher (closer to unity). The correction factor given in Eq. (25) supports this 
explanation. Lochiel and Calderbank (1964) have demonstrated that liquid drops (ethyl 
acetate drops in water and water drops in isobutanol) with larger diameter (higher Res) are 
characterized with higher correction factors. 
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Fig. 13. Effects of both superficial gas velocity uc and gas density pc (operating pressure P) on 
the correction factor f. Gas distributor: D1—perforated plate, 19 x Ø 1x10? m, triangular pitch 


The presented correction of the penetration theory is applicable in the homogeneous flow 
regime of bubble column operation. The upper boundary (so-called transitional gas 
velocity) of this flow regime was determined by the formulas of Reilly et al. (1994). 

Only in the case of two other organic liquids, carbon tetrachloride and ethylene glycol 
(reported by Oztürk et al., 1987), the presented correction method was not capable of fitting 
the experimental kia data satisfactorily. Table 1 shows that the correction method is 
applicable to organic liquids with densities pr in the range of 714-1234 kg m? and viscosities 
pL in the range of 0.327-2.94x105 Pa s. Carbon tetrachloride has a higher liquid density 
(p.71593 kg m3), whereas ethylene glycol has a higher liquid viscosity (j719.94x109 Pa s). 
In order to predict the kra values (referred to liquid volume) at high temperature conditions, 
Nedeltchev et al. (2010) used another definition of the specific gas-liquid interfacial area a 
which is referred to unit liquid volume (see Lemoine et al., 2008): 


DE (28) 
e = Ec) d, 

It is worth noting that experimental gas holdups were used for calculating those interfacial 
areas. Two different elevated temperatures T (323 and 343 K) and three different operating 
pressures P (0.1, 0.2 and 0.5) were considered, respectively. 

In the case of ellipsoidal bubbles, it is reasonable to assume that the correction factor f. is 
proportional to the interphase drag coefficient Cp and the bubble eccentricity term (E-1) 
(taking into account the distortion from the perfect spherical shape). According to Olmos et 
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al. (2003) and Krishna and van Baten (2003) the drag coefficient Cp for ellipsoidal bubbles is 
a function of Eo: 


Cp = = JEo (29) 


The drag coefficient Cp accounts for the bubble-bubble interactions and bubble-induced 
turbulence. Raymond and Zieminski (1971) reported that there is a similarity in the behavior 
of the mass transfer coefficient and drag coefficient. The decrease of the mass transfer 
coefficient corresponds to an increase of the drag coefficient (due to friction drag for small 
bubbles or form drag for large bubbles) in the case of straight chain alcohols. According to 
Raymond and Zieminski (1971) the drag coefficient is proportional to the minor semi-axis of 
the ellipsoidal bubble. The authors reported also some dependence of the mass transfer 
coefficient on the bubble shape. Wellek et al. (1966) argue that the term (E-1) can be also 
correlated to the Eötvös number (see Eq. (26)). 

Nedeltchev et al. (2010) expressed the correction factor f. for prediction of kia values 
(referred to liquid volume) as a product of both Cp and (E-1) as well as dimensionless 
temperature ratio: 


T 
f.=BCo(E-1( 2) Go 
where B is proportionality constant. 

The use of dimensionless temperature ratio is needed since the correlations for both Cp and 
(E-1) as well as Sauter-mean bubble size ds (Eq. (15)), geometrical characteristics of the 
ellipsoidal bubbles (Eqs. (10a) and (10b)) and bubble rise velocity (Eq. (5)) were derived at 
ambient temperature and thus it is not certain whether they remain valid at elevated 
temperatures. It is worth noting that Cockx et al. (1995) introduced a correction term for the 
volumetric interfacial area which is also a function of the bubble eccentricity. Talvy et al. 
(2007a) have proven that the local bubble interfacial area is sensitive to the bubble shape. For 
bubble sizes above 3x103 m they have reported correction factors above 1.1. Talvy et al. 
(2007b) have shown that the ellipsoidal shape of the bubbles seems to be significant in 
estimation of the drag coefficient. When these correlations are improved and become 
applicable at various temperatures, then the dimensionless temperature ratio in Eq. (30) 
might become useless. The substitution of Eqs. (26) and (29) into Eq. (30) yields: 


f, =0.109B Eo!” Cas) (31) 
273.15 


By means of a non-linear regressional analysis applied to 167 kra values (referred to liquid 
volume), the constant of proportionality B was estimated as 1.183. All kta values were 
successfully fitted with an average relative error of 9.3 %. 

The above-mentioned approach was also applied to slurry bubble columns operated in the 
homogeneous regime (Nedeltchev et al., 2007b). Six different liquid-solid systems 
(water/activated carbon, water/aluminium oxide, tetralin/aluminium oxide, 0.8 M sodium 
sulfate solution/kieselguhr (diatomaceous earth), ligroin (petroleum ether)/polyethylene 
and ligroin/polyvinylchloride) were considered. 85 experimental kia values were fitted 
with a mean error of 19 %. 
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5. Conclusions 


New approaches for the prediction of both gas holdups and volumetric liquid-phase mass 
transfer coefficients in gas-liquid and slurry bubble columns operated in the homogeneous 
regime have been developed. In these models the shape of the bubbles plays a major role. In 
the model for gas holdup prediction, the bubble surface and geometrical characteristics of 
the bubble (its length and height) are considered. In the model for kra prediction the same 
parameters along with the circumference of the ellipsoidal bubbles (taking part in the rate of 
surface formation) are involved. In the process of evaluation of the mass transfer coefficient, 
it is made clear that Higbie's (1935) equation should be corrected due to the non-spherical 
shape of the bubbles. New correction factors were derived. The Eótvós number Eo plays a 
major role in them. Mass transfer coefficients at elevated temperatures can also be predicted 
by this approach. 

As many as 23 organic liquids, 17 liquid mixtures and 6 liquid/solid systems were tested 
and the two major mass transfer parameters were predicted successfully. 
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7T. Nomenclature 


A = cross-sectional area of the reactor [m?] 

a = specific interfacial area (referred to dispersion volume) [m=] 
B = constant in Eqs. (30) and (31) [-] 

Cp = drag coefficient [-] 

D. = column diameter [m] 

de = bubble equivalent diameter [m] 
D; = molecular diffusivity [m2-s1] 
d, = Sauter-mean bubble diameter [m] 
e,E - bubble eccentricity (9a) or (27) [-] 

Eo = Eötvös number, Eq. (24) [-] 

fg = bubble formation frequency [s1] 

fe = correction factor [-] 

g = gravitational acceleration [m:s?] 
h - height of an ellipsoidal bubble [m] 

H - aerated liquid height [m] 
Ho = clear liquid height [m] 
kc =gas-phase mass transfer coefficient [m-s+] 
kca = volumetric gas-phase mass transfer coefficient [s1] 
kų = liquid-phase mass transfer coefficient [m:s-] 
kra = volumetric liquid-phase mass transfer coefficient [s1] 

1 = length (width) of an ellipsoidal bubble [m] 
Mo - Morton number, Eq. (14) [-] 
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Ns = number of bubbles in the dispersion [-] 
P = operating pressure [Pa] 
AP - pressure difference between the readings of both lower (at 0 m) and upper 

(at 1.2 m) pressure transducers [Pa] 
Oc = gas flow rate [m3-s-1] 
Rs = rate of surface formation [m2-s-1] 
Re = Reynolds number [-] 
Res = bubble Reynolds number, Eq. (13) [-] 
Sp = bubble surface [m2] 
Ta = Tadaki number, Eq. (12) [-] 
T -operating temperature [K] 
te = contact time [s] 
up = bubble rise velocity [m:s1] 
uc = superficial gas velocity [m-s+] 
Utrans = transition gas velocity [m-s-1] 
Vs = bubble volume [m3] 
Vc = gas volume [m3] 
Vi = liquid volume [m3] 
Viotal = dispersion volume [m3] 
We = Weber number [-] 


8. Greek letters 


£c = gas holdup [-] 

pL = liquid viscosity [Pa-s] 
Pc= gas density [kg-m3] 
pctef = reference gas density (air at ambient conditions) [kg:m?] 
pL = liquid density [kg:m:?] 
OL = liquid surface tension [N-m-1] 
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